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I , I The PREX experiment at Jefferson Laboratory measures the neutron radius of 

208Pb with parity violating electron scattering in a way that is free from most 
I strong interaction uncertainties. Tlie 208Pb radius lias important implications 

for neutron rich matter and the structure of neutron stars. We present first 
PREX results, describe future plans, and discuss a follow on measurement of 
the neutron radius of 48Ca. We review radio and X-ray observations of neutron 
star masses and radii. These constrain the equation of state (pressure versus 
density) of neutron rich matter. We present a new energy functional that is 
^— H simultaneously fit to both nucleaj: and neutron star properties. In this approach, 

^ neutron star masses and radii constrain the energy of neutron matter. This 

avoids having to rely on model dependent microscopic calculations of neutron 
l/^ matter. The functional is then used to predict the location of the drip lines 

and the properties of very neutron rich heavy nuclei. 

cn 
o 

Compress matter to high densities (10^^+ g/cm^) and electrons are cap- 
. . turcd to make neutron rich matter. This material is at the heart of many 

_ ^ fundamental questions in nuclear physics and astrophysics. 



>< 



1. Introduction 



• What are the high density phases of QCD? 

• Where did the chemical elements come from? 

• What is the structure of many compact and energetic objects in 
the heavens, and what determines their electromagnetic, neutrino/ 

and gravitational-wave^ radiations? 

We are in Interested in neutron rich matter over a great range of density 
and temperature were it can be a gas, liquid, solid, plasma, liquid crystal 
(nuclear pasta), superconductor {T^ = 10^^ K), superfluid, color supercon- 
ductor. The liquid crystal phases are known as nuclear pasta.''''* Pasta is 
expected at the base of the crust in a neutron star and can involve complex 
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shapes such as long rods ("spaghetti") or flat plates ( "lasagna" ) , see Fig. 
1. Pasta may scatter neutrinos in core collapse supernovae because it has 
sizes comparable to neutrino wave lengths.^ 




Fig. 1. Surfaces of proton density for a pasta configuration of neutron rich matter at a 
baryon density of 0.05 fm~^. This is from a scmiclassical molecular dynamics simulation 
with 100,000 nucleons.*^ 

In Section 2 we describe a precision laboratory experiment called PREX 
to measure the neutron radius of ^"^^Pb. In Section 3 we describe a new 
energy density functional that is fit to properties of finite nuclei and the 
maximum mass and radius of neutron stars. We conclude in Section 4. 

2. The Lead Radius Experiment (PREX) 

Neutron rich matter can be studied in the laboratory. We give one example, 
the Lead Radius Experiment (PREX)^ accurately measures the neutron ra- 
dius in ^"^^Pb with parity violating electron scattering.* This has important 
implications for nuclear structure, astrophysics, and atomic parity violation. 
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Parity violating electron scattering provides a model independent probe 
of neutron densities that is free from most strong interaction uncertainties. 
This is because the weak charge of a neutron is much larger than that of a 
proton.^ Therefore the boson, that carries the weak force, couples pri- 
marily to neutrons. In Born approximation, the parity violating asymmetry 
Apy, the fractional difference in cross sections for positive and negative he- 
licity electrons, is proportional to the weak form factor. This is the Fourier 
transform of the weak charge density pw{r), see Fig. 2, and is closely re- 
lated to the neutron density. Therefore the neutron density can be extracted 
from an electro- weak measurement.^ Note that the Born approximation is 
not valid for a heavy nucleus and coulomb distortion effects must be in- 
cluded. However, these have been accurately calculated. Many details of 
a practical parity violating experiment to measure neutron densities have 
been discussed in a long paper. ^ 

The Lead Radius Experiment (PREX) at Jefferson Laboratory^ mea- 
sures the parity violating asymmetry Apy for 1.05 GeV electrons elastically 
scattered from ^"^Pb at laboratory angles near five degrees. The first mea- 
surement yielded Apy — 0.656 ± 0.060 (statistical) ±0.014 (systematic) 
ppm.'' From this the weak charge density has been deduced, sec Fig. 2, and 
the rms neutron radius i?„ minus proton radius Rp for ^°*Pb was found to 
be Rn — Rp = 0.33^^018 f™- ref.^^ for details of this analysis. A second 
PREX run is now approved to accumulate more statistics and reach the 
original goal of determining i?„ to 1% (±0.05 fm). 

In addition to PREX, many other parity violating measurements of 
neutron densities are possible, see for example. Measuring ii„ in ^^Ca is 
particularly attractive. With only 48 nuclcons, microscopic coupled cluster 
calculations,^"^ or no core shell model calculations,^'* are feasible for ''^Ca 
that are presently not feasible for ^°^Pb. Note that these microscopic cal- 
culations have important contributions from three nucleon forces. This will 
allow one to make microscopic predictions for the neutron density and relate 
a measured i?„ to three nucleon forces. 

The neutron radius of ^°^Pb, i?„, has implications for astrophysics. 
There is a strong correlation between i?„ and the pressure of neutron mat- 
ter P at densities near 0.1 fm""^ (about 2/3 of nuclear density). A larger 
P will push neutrons out against surface tension and increase Rn- There- 
fore measuring i?„ constrains the equation of state (EOS) — pressure as a 
function of density — of neutron matter. 

Recently Hebeler et al.*^ used chiral perturbation theory to calculate 
the EOS of neutron matter including important contributions from very 
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Fig. 2. Weak charge density pw{^) of ^"^Pb that is consistent with the PREX result 
(solid black line).^^ The brown error band shows the incoherent sum of experimental 
and model errors. The red dashed curve is the experimental (electromagnetic) charge 
density p^^. 



interesting three neutron forces. We have some information on isospin 1/2 
three nucleon forces from mass 3 nuclei ("^He, '^H) and proton-deuteron 
scattering. However, our experimental information on three neutron forces 
is limited. From their EOS, they predict Rn — Rp = 0.17±0.03 fm. Here Rp 
is the known proton radius of ^°®Pb. Monte Carlo calculations by Carlson et 
al. also find sensitivity to three neutron forces. "'^^ Therefore, measuring i?„ 
provides an important check of fundamental neutron matter calculations, 
and constrains three neutron forces. 

The EOS of neutron matter is closely related to the symmetry energy 
S. This describes how the energy of nuclear matter rises as one goes away 
from equal numbers of neutrons and protons. There is a strong correlation 
between f?„ and the density dependence of the symmetry energy dS/dn, 
with n the baryon density. The symmetry energy can be probed in heavy ion 
collisions. "'^^ For example, dS/dn has been extracted from isospin diffusion 
data"'^^ using a transport model. 

The symmetry energy S helps determine the composition of a neutron 
star. A large 5, at high density, implies a large proton fraction Yp that 
will allow the direct URCA process of rapid neutrino cooling. If i?„ — R^ 
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is large, it is likely that massive neutron stars will cool quickly by direct 
URCA.^° In addition, the transition density from solid neutron star crust 
to the liquid interior is strongly correlated with i?„ — Rp?^ 

Finally, the correlation between _R„ and the radius of a neutron star rjvs, 
see for example Ref.,^^'^^ is also very interesting.^'* In general, a larger i?„ 
implies a stiffer EOS, with a larger pressure, that will also suggest rjvs is 
larger. Note that this correlation is between objects that differ in size by 
18 orders of magnitude from i?„ w 5.5 fm to rjv5 ~ 10 km. 




I ^ ^ "'1 ■ =^ 

6 8 10 12 14 16 
Radius R (km) 

Fig. 3. Neutron star mass- radius relation of SLy4,25 UNEDFO^" and SV -min.^'^ The 
uncertainty band for SV-min is shown. This band is estimated by calculating the covari- 
ance ellipsoid for the mass M and the radius R at each point of the SV-min curve as 
indicated by the ellipsoid. Also depicted (dotted lines) are uncertainty limits for TOV- 
min. 



3. Energy density functional for nuclei and neutron stars 

Energy density functionals (EDF) are often fit to properties of finite nuclei. 
However, this can leave many properties of bulk neutron matter poorly 
constrained. For example the SV-min functional has large uncertainties in 
its predictions of neutron star mass versus radius, see Fig. 3. Alternative 
optimization protocols are based on considering both experimental data on 
finite nuclei and theoretical pseudo-data on nuclear and neutron matter. 
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For example, the SLy4^^ EDF was fit to both nuclei and the variational 
Fermi hypernetted chain calculations of neutron matter. Chamel et al. 
calculated the maximum mass of neutron stars for three EOS fit to finite 
nuclei and microscopic neutron matter calculations,^^ and such analysis 
helped them to so select one EDF (BSk21) that does well on both kinds of 
data. 

However, it is to be noted that current microscopic calculations of neu- 
tron matter with phenomenological interactions are highly model depen- 
dent; in particular, the impact of the poorly known three-neutron forces 
(especially their T = 3/2 component) seems to be crucial. ^^'^^ In this re- 
spect, calculations based on the chiral cfToctive field theory with all many- 
neutron forces are predicted to n3lq30,31,36 provide useful guidance, 
at least at sub-saturation densities. At high density, however, four or more- 
nucleon forces^^'^® may also be important. Indeed, at this time, our ability 
to calculate neutron matter properties at high densities is fundamentally 
limited. 

Recently we fit a new energy density functional, TOV-min, to both 
properties of finite nuclei and to the maximum mass of a neutron star M^ax 
and the radius of a 1.4 Mq solar mass neutron star"^^ (see also rcf.'^^). This 
functional makes much more accurate predictions for neutron star mass 
versus radius, see Fig. 3. The functional, with neutron matter properties 
constrained by neutron star data, then may make more accurate predictions 
for very neutron rich nuclei. For example. Fig. 4 shows the location of the 
drip line for a large range of nuclei. 



4. Conclusions: neutron rich nuclei and neutron stars 

Neutron rich matter is at the heart of many fundamental questions in Nu- 
clear Physics and Astrophysics. What are the high density phases of QCD? 

Where did the chemical elements come from? What is the structure of many 
compact and energetic objects in the heavens, and what determines their 
electromagnetic, neutrino, and gravitational-wave radiations? 

We described the Lead Radius Experiment (PREX) that uses parity 
violating electron scattering to measure the neutron radius in ^°*Pb. This 
has important implications for neutron stars and their crusts. We also de- 
scribed a new energy functional fit to the properties of finite nuclei and to 
the masses and radii of neutron stars. This functional then predicts prop- 
erties for very neutron rich nuclei. 
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Fig. 4. Two-nucleon drip line of TOV-min (dark blue)''^ compared to the uncertainty 
band from Ref.^* (2012 Benchmark) obtained by averaging results of six different EDF: 
SkM*, SkP, SLy4, SV-min, UNEDFO, and UNEDFl. The dashed grey grid lines show 
the magic numbers known around the valley of stability (20, 28, 50, 82, 126) as well as 
the predicted regions of enhanced shell-stability in superheavy nuclei around N = 184 
and 258. All 767 experimentally known even-even isotopes are shown with the stable 
nuclei indicated as black squares and the radioactive nuclei as green squares. 
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